To appear in Astrophys J. March 10, 2002 

Ambipolar Drift in a Turbulent Medium 

Ellen G. Zweibel 

JILA & Department of Astrophysical & Planetary Sciences, University of Colorado, Boulder CO 

^ 80309 

O 

^ ■ ABSTRACT 

> 

o _ 

^^ ' The interstellar magnetic field strength and density are observed to be correlated, 

\f~^ • but there is a large dispersion in this relation. In particular, the magnetic field is often 

observed to be weaker than expected. At low ionization fraction, ion-neutral drift, or 
^S| ■ ambipolar diffusion, permits slip of the field relative to the neutral gas and tends to 

J^ , make the field strength more uniform, but is thought to be too slow to explain the obser- 

ve I vations. The purpose of this paper is to show that ion-neutral drift is significantly faster 

!^ I in a turbulent medium than in a quiescent one. We suggest that this fast ambipolar 

f^ ' diffusion can explain the surprisingly low magnetic fieldstrengths sometimes observed 

in dense interstellar gas. 
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52 I 1. Introduction 

The interstellar magnetic fieldstrength and gas density are observed to be correlated (Troland 
rS I Sz Heiles 1986, Crutcher 1999). This relationship is thought to arise from so-called ideal magne- 

C^ I tohydrodynamic (MHD) processes, in which the field is frozen to the ambient medium. Under 

ideal MHD conditions, if the ratio of mass to magnetic flux were everywhere constant, the slope 
q = dlogB/dlogp of the fieldstrength - density correlation would be unity for compression normal 
to B, zero for compression parallel to B, and 2/3 for isotropic compression. The observed value of 
q is approximately 0.5, which is consistent with equilibrium models of self gravitating clouds which 
evolved under conditions of frozen flux (Mouschovias 1976, Tomisaka, Ikeuchi, & Nakamura 1988). 

However, a number of observations and upper limits on magnetic fieldstrength suggest that 
B oc p^'^ is more an upper envelope than a scaling law. This is true both in atomic and molecular 
gas (Bourke et al. 2001, Crutcher 1999, Heiles 2001a, Heiles & Troland 2001) Although the number 
of measurements is small, and the field may be underresolved in some cases (Brogan &: Troland 
2001), the trend towards weak fields is clear. A different line of argument comes from numerical 
simulations of molecular clouds: Padoan & Nordlund (1999) claim that models with weak fields 
replicate the observations better than models with strong fields. 



Weak fields are difficult to reconcile with ideal MHD. This is particularly so in the case of 
turbulent molecular clouds. Giant molecular clouds (GMCs) are about 40 times more dense than 
the mean interstellar gas, suggesting that the magnetic field should be 6-7 times stronger than the 
mean field. If GMCs are formed by flows parallel to the field then of course the field would not 
be strengthened, but the mean field is not strong enough to resist compression and collimate the 
flows. Furthermore, Mestel (1985) has pointed out that accumulation of the mass of a GMC by 
compression in one dimension requires organized motion of the gas over nearly one kiloparsec, the 
origin of which is difficult to understand. 

Thus, we seek an explanation for the weakness of the magnetic fieldstrength - gas density 
correlation beyond the scope of ideal MHD. At the very largest lengthscales, ideal MHD should be 
an excellent approximation. Moving downwards in scale, the first non-ideal effect encountered, at 
least at low ionization fraction, is ion-neutral drift, or ambipolar diffusion. At the ambipolar scale, 
which is many orders of magnitude larger than the resistive scale, the magnetic field and plasma 
become decoupled from the neutral material. This makes it possible to change the mass to flux 
ratio without altering the magnetic topology. Ambipolar drift has been invoked as the primary 
magnetic flux transport mechanism in dense, star forming gas since the classic paper by Mestel &: 
Spitzer (1956). However, it is thought to be too slow to be an effective transport mechanism in 
diffuse gas (see §2 for quantitative estimates). 

Interstellar gas is turbulent. Turbulent diffusion of quantities such as heat and angular mo- 
mentum is often invoked in astrophysics as a mechanism for enhancing transport rates above their 
kinetic theory values, which are usually very slow. Turbulence enhances diffusion rates by mix- 
ing the relevant quantity to the small scales at which molecular diffusion operates. This leads to 
a mixing time which is approximately the eddy turnover time, and is nearly independent of the 
molecular diffusivity. 

Whether turbulence enhances the resistive decay rate of a magnetic field is unclear, because 
there is substantial evidence that magnetic forces resist stretching the field sufficiently to mix it to 
the tiny scales at which resistivity operates (Cattaneo & Vainshtein 1991, Cattaneo 1994, Cattaneo, 
Hughes, & Kim 1996). This paper addresses a different question, namely, whether turbulence in a 
weakly ionized gas can transport magnetic field with respect to the neutral matter, without resistive 
dissipation necessarily coming into play. Because the ambipolar drift scale is much larger than the 
resistive scale, the feedback effects which can suppress turbulent resistivity are far less dramatic, 
although they cannot always be ignored. 

We use analytical methods to calculate the mixing rate. Numerical study of mixing requires 
that numerical diffusion of both field and fluid be very well controlled. Analytical calculations are 
useful for initial exploration of some of the basic mechanisms. 

In §2, we introduce the physical model, derive an equation for the evolution of the mass to 
magnetic flux ratio in a weakly ionized medium, estimate relevant timescales, and establish an 
initial condition. In §3, we estimate the turbulent diffusion rate based on mixing length theory. 



quantify this estimate with an exact calculation of transport by a wave, and argue that enhanced 
diffusion requires the introduction of small scales as well as bulk advection. In §4 we develop a 
model based on exponential stretching and shrinking. This model leads to a flux redistribution 
rate which is comparable to the eddy turnover rate, as expected in turbulent diffusion problems. 
However, the model is two dimensional, and the fluid motions are prescribed without specifically 
allowing for magnetic forces. In §5, we consider magnetic feedback. In §6, we apply the model to the 
interstellar medium, and discuss the astrophysical constraints imposed by the dynamics. Section 
7 is a summary and conclusion. Sections 3 and 4 are relevant to general mixing problems, such as 
turbulent diffusion of a passive scalar, and the reader who is mainly interested in the astrophysical 
implications could go directly to §6. 



2. Basic Equations and Model 

2.1. Equations for Ambipolar Drift 

We consider a weakly ionized medium with magnetic field B, mass density p, and ion mass 
density pi <^ p. We are interested in timescales much longer than the ion-neutral collision time Tin, 
in which case the ion-neutral drift vd = v; — Vn is well approximated by 



(V X B) xB 

VD = -. Tin, (1) 

Anpi 



(Shu 1983). 



The magnetic field evolves according to the magnetic induction equation for a perfectly con- 
ducting medium 

(9B 

— = V X (V; X B), (2) 

Replacing Vj with vd + Vn, and approximating Vn by the center of mass velocity v, we rewrite the 
induction equation as 

= V X (v X B) + V X (vD X B), (3) 

ot 

where vd is given by eqn. (1). The first term on the RHS of eqn. (3) can be expanded using the 

identity V x (v x B) = B • Vv — v • VB — BV • v. Then, using the continuity equation 

^ = -vVp-pV-v (4) 

we derive an evolution equation for the magnetic field to density ratio B/p 

d B B B 1 , , 

hv- V— = — • Vv-F-V X (vD X B). (5) 

dt p P P P 

The left hand side of eqn. (5) is the comoving, or convective, time derivative of B/p. The first term 
on the right hand side represents stretching of the fieldlines, and is a consequence of the frozen 



flux condition. The second term on the right hand side represents the evolution of B/p caused by 
ambipolar drift. 

We now restrict ourselves to two dimensional, incompressible flows perpendicular to a straight 
magnetic field. This geometry captures the main effects under study, and is consistent with the 
nature of turbulence in a strong, well ordered magnetic field (Strauss 1976, Sridhar & Goldreich 
1994, Goldreich & Sridhar 1997). We recognize that interstellar turbulence is frequently observed 
to be highly supersonic, and thus cannot be entirely incompressible. We expect that compressible 
turbulence to result in magnetic flux transport just as incompressible turbulence does, but that 
magnetic feedback on the turbulence is stronger in the compressible case. 

Under the assumptions of incompressibility and two dimensionality, the line stretching term 
vanishes, and eqn. (5) reduces to 

d B B I B'^ 

+ v.V- = -V-^T„,V5, 6 

ot p P P 47rp 

where B now represents the amplitude of the magnetic field, Tni = nnpn/pi is the neutral-ion 

collision time, and we have used eqn. (1). Equation (6) is close to an advection - diffusion equation 

for the flux to mass ratio Q = B/p 

^ + vVQ = iv-AVB, (7) 

ot p 

provided that we define the diffusivity A as Tniv\, where B/{4ttp)^'^ is the Alfven speed. 

The diffusion of Q is nonlinear in the sense that A = X{Q,p). This nonlinearity can produce 
sharp fronts along surfaces where B vanishes or is tightly folded (Brandenburg &: Zweibel 1994), 
similar to fronts created by nonlinear thermal conduction (Zel'dovich & Raizer 1966). Resistive 
diffusion in these current sheets alters the mass to flux ratio as well as changing the magnetic 
topology (Brandenburg & Zweibel 1995, Zweibel & Brandenburg 1997). In this paper we assume 
that the relative variation of B/p is so weak that nonlinear effects play only a minor role in 
ambipolar drift, and B remains nonsingular. 

Equation (7) can be used to derive an equation for the rate of change of B within a comoving 
volume V of fluid; that is, a fluid element of fixed mass. We have 

— / d^xpQ = / d^x—pQ + / d^xpQv ■ n, (8) 

dt Jy Jv Ot Js 

where S is the surface of V. The first term on the RHS of eqn. (8) is the Eulerian change of B and 
the second term accounts for the motion of V. Using eqns. (4) and (7) to expand the first term on 
the RHS and applying Gauss' theorem yields 



, , d-^xpQ= / d^xXh-VB. (9) 

at Jv Js 

Since the mass within V is constant, eqn. (9) shows that the flux to mass ratio within a volume 

moving with the center of mass velocity decreases if the magnetic field decreases outward on the 

surface of the volume. 



2.2. Timescales 

In the absence of flow, the characteristic diffusion time for a magnetic field of representative 
strength B and scale length L ~ \B/VB\ is 

tdo = ^ = ^ = — ,• (10) 

^ V j^Tni Tni 

where ta = L/va is the global Alfven time. 

Expressing t^o in physical units reveals the magnitude of the timescale problem. We take 
Tni from Draine, Dalgarno, & Roberge (1983); when the ratio of ion to neutral atomic weight 
Ai/An > 1, Tni = 6.7 X W^n~^ s. The Alfven speed va = 2.2 x W^B^/{nnAn)^/'^ cm s~^. The 
diffusivity A is then 3.2 x 10^^i3^/(yl„n^Xj) cm^ s~"^, where Xj = ni/un is the ionization fraction. 
The drift time is 

t.0 = 3.1xl020^%^s, (11) 

where A^20 = finL is the column density in units of 10^*^ cm^^. For example, the systems reported 
by Heiles (2001) have A'^20 a few tenths to a few, with B^ typically 3. If we take An = 1.4 for gas 
of cosmic composition and Xi = 10~^, which is probably a conservative lower limit, then we find 
trfo is of order 10^ yr for these systems. Since this is much more than the 10^ - 10^ yr expected 
lifetime of an interstellar cloud, the diffusion rate must be enhanced by a factor of 10 - 100 in order 
to explain the flatness of the B — n relation. 

Mixing by eddies in the neutral gas requires that the magnetic field be frozen to the flow. The 
degree of freezing is measured by the ambipolar Reynolds number Rad, which is large under frozen 
in conditions (Zweibel & Brandenburg 1997). For eddies of characteristic size I and speed vt, 

Rad ^ ^. (12) 

Thus, the field is frozen to the turbulent fiow for 

j) {^) > ^. (13) 

LJ \vaJ ta 

Since t^iIta is expected to be small, while vt/vA is order unity, eqn. (13) implies a fair degree of 
dynamic range for the sizes of eddies that can mix the field. Equation (12) can also be written as 

Rad = 4 — ' (14) 

^A ''""* 

which shows that if the fiow is at or above equipartition with the field {vt > va), the field is frozen 
in if the neutral-ion collision time is less than the eddy turnover time. 



2.3. Initial Condition 

It will be useful in the following analysis to have a definite model for B. We take as an initial 
condition 

B{x,y,0) = Boo + ^BQX^, (15) 

where Bqq and Bq < are constants. We will define A using Bqq for Bq in quantitative examples. 
We view eqn. (15) as the first two terms in a Taylor expansion of a magnetic field which peaks at 
X = 0, and will assume x/L ^ 1. 

We first consider pure diffusion. Motivated by the initial condition (15), we seek solutions of 
eqn. (7), with v = 0, of the form 

B{x,y,t) = Bo{t) + ^B'y. (16) 

Substituting eqn. (16) into eqn. (7) and using eqn. (15), we find 

Boit) = Boo + Bo\t. (17) 

Equation (17) predicts that the peak field decreases by a factor of two on a timescale 

tdo = -Soo/(2ASo). (18) 

Equation (18) agrees with eqn. (10) if we define the magnetic lengthscale L by 

3. Diffusion in the Presence of Waves 

We begin with a mixing length argument. Consider a magnetic flux tube of width a which is 
carried by a random flow u a distance / down the gradient of B. The field in the tube diffuses into 
the ambient medium; thus the motion causes net transport of B. The transport is most effective 
when the diffusion time across the tube is comparable to the advection time 

(20) 

because if a^/A ^ l/u, the flux tubes return to their original positions with nearly the same value 
of the field, while if a^/A <^ l/u the field diffuses too quickly to be advected by the flow. 

Advection spreads B over a distance / in a time l/u. In this time, B would spread diffusively 
over a distance (A//u)^". By eqn. (20), this distance is just a. Advective mixing accelerates the 
transport of B only if a/l < 1, meaning that the motion consists of thin fingers which travel much 
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T 




u 



further than their widths (see Ottino 1989 for discussion of mixing by fingers, or tendrils). Such 
fingers are not seen in models of Alfvenic turbulence, and it is not clear that they would form in 
a weakly stratified gas such as the interstellar medium. Interpreted more broadly, the argument 
presented here shows that turbulent mixing requires more than just advection and dispersal; it also 
requires the formation of small scales. 

Now, we quantify this result. Weak turbulence theory, in which turbulence is modelled as a 
superposition of randomly phased waves, can be used to compute the rate of turbulent diffusion (eg. 
Moffatt 1978, Gruzinov & Diamond 1994). In this so-called quasilinear approach, one partitions 
quantities into mean and fluctuating parts and calculates the average effect of the fluctuations on 
the mean part. We used this method in a previous study of ambipolar diffusion (Zweibel 1988). 
In the present problem it is possible, as well as instructive, to solve the induction equation exactly 
instead of averaging it. This confirms the argument given in §2. 

We introduce a periodic flow in the x direction 

V = xusinLJtsinfcy. (21) 

Motivated by the initial condition eqn. (15), we try a solution of the advection-diffusion equation 
(7) of the form 

B{x,y,t) = Bo{t) + -B'o + Bi{t)x sin ky + B2{t) cos 2ky. (22) 

Substituting eqn. (22) into eqn. (7), using eqn. (21), and equating like powers of x and Fourier 
harmonics of y leads to a set of coupled ODEs for the functions Bq, Bi, and i?2. The solution for 
Bq, which follows the decay of the peak magnetic field, is 



^0 = ^00 + ^Bq t 



n" 2 



2 (tj2 + r2) 



,2 



r / sin2ci;t\ sin^ ut 



UJ 



2-— ^2(1-6 COSWt) 



uT 



U) 



2 + r2 



—Vf 

e sm ut 



(23) 



where F = Xk'^. 



The maximum decay rate occurs when the motion given by eqn. (21) is coherent over many 
wave periods. The long time behavior of Bq is then given by 

Bo = Boo + (A + Ai) B'ot, (24) 

where 

Fii2 

^* ^ 4(c.2 + F2) (2^) 



represents diffusion brought about by advective transport. Equation (25) closely resembles the 
turbulent diffusivity calculated from quasilinear theory (Moffatt 1978). Maximizing A^ with respect 
to k, we find that the maximum occurs for oj = F, as we asserted in the mixing length argument 
following eqn. (10), and is 

^t,max = g^^2' ^ ) 

where we have replaced F by A/c^. If we express u in terms of the maximum fluid displacement 
a = ujijj and take the ratio of Xt,max to A, the result is 

A " 8 ■ ^^' 

Equation (27) shows that the diffusion rate is appreciably enhanced by waves only if fca ^ 1, 
meaning that the flow consists of long, thin streamers (see also Press &: Rybicki 1981). We reached 
the same conclusion from mixing length theory. The missing ingredient is stretching and shrinking 
of scales, an intrinsic feature of turbulent flows which we incorporate in the next section. 



4. Stagnation Point Flow 

Hyperbolic stagnation point flow is a particularly tractable example of a flow with exponen- 
tial shrinking and stretching. At hyperbolic stagnation points, the fluid flow converges in one (or 
two) directions and diverges in the other direction(s), while maintaining incompressibility. It is 
well known that diffusion is accelerated in the vicinity of stagnation points, due to the shrinking 
of scales in the convergent directions (Moffatt 1978, Zweibel 1998), while Zel'dovich et al. (1984) 
demonstrated dynamo action by a random ensemble of stagnation points. The role of hyper- 
bolic stagnation points in the mixing of scalar fields in turbulent flows has recently been reviewed 
(Shraiman & Siggia 2000) with emphasis on the development of intermittency, and the high order 
moments of the distribution of concentrations. 

The advection-diffusion problem for fields of the form (15) is exactly soluble for stagnation 
point flow. In subsection (4.1) we compute the effect of a single stagnation point. In subsection 
(4.2) we superimpose the effects of a random ensemble of stagnation points, and in subsection (4.3) 
we discuss the relationship between the stagnation point model and turbulent flow. Our model is 
not intended to be a full theory of turbulence, but merely illustrative. 



4.1. A Single Stagnation Point 

We consider two dimensional, incompressible flow near a stagnation point at (ax, ay). For fields 
of the form given by eqn. (15), we require a^/L <^ 1, where L is given by eqn. (19). The flow is 

Va: = --f {x - a^:) , (28) 

iJy=l{y- ay) , (29) 



where 7 is a constant. It is straightforward to integrate equations (28) to find the position at time 
t of a fluid parcel which is at position (xq, yo) at t = 

X = Qx + (xo - ax) e"^* (30) 

y = ay + {yo - ay) e<K (31) 

The initial coordinates in terms of the coordinates at time t are 

Xq = ax + {x - ax)e'* (32) 

Vg = ay + [y - ay)e-~'^^ ■ (33) 

We now compute the effect of this stagnation point flow on the diffusion of the magnetic field. 
With eqn. (15) as the initial condition, we look for a solution of the form 

B{x,y,t) = Bo{t) + ]^Blxl (34) 

where Xo{x,y,t) is given by eqn. (32). Substituting eqn. (34) into the advection-diffusion equation 
(7) yields 

^0 = XB^e^^K (35) 

B is of the form (34) exist for two reasons. First, xq is a constant of the motion, so for any function 
f{xo) 

^^ + v.V/(xo) = 0, (36) 

where v is given by eqns. (28). Second, eqn. (32) shows that xq is a linear function of x and y, so 
V^xo is only a function of time, and is independent of the stagnation point location {ax, ay). 

The solution of eqn. (35) which fits the initial conditions is 

Bo{t)=B,o + ^Bl{e^'^'-l). (37) 

The location of the peak field evolves in time to ax{i — e^^^), but this is irrelevant because the 
density p is shifted by the same amount. It is only diffusion which affects the mass to flux ratio. 

According to eqn. (35), the peak field has dropped to half its value in the time t.y 

t^ = -Lin{l + 2jtdo), (38) 

27 

where tdo is the diffusion time in the static case, defined in eqn. (18). Equation (38) shows that 
the diffusion time depends only logarithmically on the diffusivity A. This arises because of the 
exponential growth of the magnetic field gradient, as seen in eqn. (35). 

The quantity 2^tiiQ which appears in the logarithm in eqn. (38) is, however, a large number. 
If we identify 7"-*^ ~ l/vt with an eddy turnover time r^, and use eqns. (10) and (12), we can write 

2jtdo ~ 2 —^ ^2Radj2- (39) 



10 



At the time t^, the x component of the drift velocity vdx is half the flow velocity Vx = 'jx. 
Beyond this time, the magnetic field is not well coupled to the flow. 



4.2. An Ensemble of Stagnation Points 

We compute the evolution of the magnetic field carried by Lagrangian fluid elements under 
sequences of stagnation point flows oriented in random directions, each one of which endures for a 
time T (see Childress &: Gilbert 1995 for a general discussion of these so-called renewing flows). 

We take the flow during time (n — l)r < t < nr to be 

Vnx = IfJ'uX + 7(1- fJ'l) y, (40) 

1 /9 

Vny = 7 (1 - ^n) ^ " It^nV, (41) 

where — 1 < ^n ^ 1 and n > 1, and for simplicity we have taken all flows to have the same strength 
7. Equation (40) reduces to eqn. (28) if ^ = —1. These flows are curl free, and hence not of the 
most general possible type. However, in the neighborhoods of stagnation points, vorticity leads to 
changes of scale at an algebraic rather than exponential rate, complicating the mathematics while 
having little effect on the diffusion rate (Zel'dovich et al. 1984, Zweibel 1998). Thus, vorticity is of 
secondary importance to our problem, and we omit it here, although interstellar turbulence almost 
certainly possesses vorticity. 

Let the {x,y) coordinates of a fluid parcel at time (n — 1)t be rn_i. Then at time nr, the 
coordinates rn can be written as 

rn = AnTn_i, (42) 



where the matrix An is 



cosh 7r + ^„ sinh 7r (l — /U^) sinh7r 



y (l — /U^) sinh7r cosh 7r — ^„ sinh 7r J 
Inverting eqn. (42) yields rn_i in terms of rn 

rn-l = An^Tn. (43) 

Successive backwards iteration of eqn. (43) yields the initial position ro in terms of the coordinate 

TN = r(iVT,ro) 

ro = Ai-i-A2-^..AN-^-rN. (44) 

At times nr < t < (n + 1)t, the position r(t) is related to the postion at time nr by an equation 
similar to eqn. (42), where in the matrix A we replace ^t by j{t — nr). 

Since the matrices An are independent of the spatial coordinates, ro is a linear function of 
FN, or, more generally, r(t). This means that V^Xg is a function only of time. It follows that the 
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solution B of the advection- diffusion equation (7) can still be written in the form of eqn. (34), and 
that the rate of diffusion increases with time at the same rate as V^Xq. 

In order to estimate the rate of increase of V^Xq, we generated random sequences of /i„ and 
calculated ro from rjvj according to eqn. (44). We used the result to calculate V^Xq as a function of 
n, or equivalently, as a function of time, since n corresponds to the time nr. The only adjustable 
parameter in these calculations is 'jt, which measures the coherence of the flow, here the renewal 
interval in units of the stretching rate. We expect ^t to be 0{1). 

Figure (1) shows the dimensionless diffusion rate as a function of iteration number in the case 
7r = 0.5. On average, the increase in diffusion rate is well fit by an exponential, and is about 84% 



15 - 




Fig. 1. — The natural log of the increase in diffusion rate, or amplification factor, with time, or the 
number of iterations. Each solid curve is the average of 50 independent realizations of the iteration 
process. The dashed line, which has slope 0.838, is the exponential with the same final value as 
the average of the curves. The maximum amplification factor possible would occur if all stagnation 
points had inflow along the x axis, and would have a slope of 1 in these units. 

the rate of increase for coherent stagnation point flow given in eqn. (28). After 20 renewals the 
average diffusion rate is more than 10^ times larger than its original value. This is much larger than 
the 2-3 orders of magnitude that we estimated in §2.2 as required to explain the B — n relation. 
The mean amplification rate is relatively insensitive to the coherence parameter jt, being about 
75% of maximum if 'jt = 0.1 and about 87% of maximum if jt = 1. 

However, there is substantial dispersion about the mean. Each solid curve represents 50 
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independent sequences of iterations, and Fig. (1) shows differences between them. The standard 
deviation within each set of 50 sequences is typically about 40% of the mean, and the amplification 
factors for single sequences rarely grow exponentially. This is illustrated in Figure (2). The large 




Fig. 2. — The natural log of the increase in diffusion rate, or amplification factor, with time, or the 
number of iterations. Each solid curve is a single independent realization of the iteration process. 
The curves, which correspond to the first 5 members of a larger ensemble, show the intrinsic 
variability of the amplification process. 

standard deviation suggests that the diffusion rate in this model is highly intermittent, which 
is characteristic of turbulent mixing. Additional evidence of intermittency is seen in the PDFs 
of the distribution of amplification factors, shown in Figure (3) for two different values of ^t. 
The maximum possible amplification rate of 20 imposes a cutoff on the high amplification side of 
the curve for ^t = 1.0; the PDF for the case ^t = 0.1 is more symmetrical because the mean 
amplification rate is well below the maximum. 



4.3. The Stagnation Point Model and Turbulent Flow 

The hyperbolic stagnation point model achieves fast diffusion by increasing the magnetic field 
gradient at, on average, nearly exponential rates. Chaotic fiow achieves fast stretching and shrinking 
without hyperbolic stagnation points because the trajectories of neighboring fiuid particles separate 
at an exponential rate. Small differences in the rates of exponential change lead to highly intermit- 
tent distributions of scalar quantities, which we saw refiected in the stagnation point model through 
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Fig. 3. — Normalized frequency distributions, or PDFs, of the natural logarithms of the amplifi- 
cation factors for 10 iterations at ^t = 1.0 (solid curve) and 100 iterations at 7r = 0.1 (dashed 
curve), so that the curves correspond to the same total time. Each curve is based on 10^ random 
sequences of iterations. 



the wide dispersion of amplification rates (see Figures 2 & 3). The maximum rate of stretching of 
a fluid element labelled by its initial position xq is given by the Lyapunov exponent A(xo) 



A (xq j = max lim sup — In | — — • Gq | , 
eo i^oo t axo 

where the eg are the set of all possible unit vectors (see Childress &: Gilbert 1995). 



(45) 



In order to make it plausible that the diffusion rate is enhanced by exponential shrinking and 
stretching in a flow, we imagine that the diffusivity A is so small that we can ignore it. In this 
limit, the solution of the advection equation for the initial condition (15) is 



B (x, t) = B (xo, 0) = Soo + ^B'^xl (x) . 



(46) 



If we restore diffusion, it appears in the advection-diffusion equation (7) as the term 

XB'q (xqV^o + Vxo • Vxo) . (47) 



XB'JSJ^^ 



The quantity Vxq • Vxq on the right hand side of eqn. (47) is the square of the inverse of the 
stretching rate. This suggests heuristically that the diffusion rate grows exponentially. 
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Examples of 2D, chaotic, spatially periodic flow with a single lengthscale, including maps of 
Lyapunov exponents and other measures of chaos, are given by Galloway &: Proctor (1992), Ponty 
et al. (1993) and Cattaneo et al. (1995). 



5. Dynamical Feedback with Extension to 3D 

In §4, we prescribed a strictly 2D flow, and neglected feedback by Lorentz forces. In fact, 
by virtue of eqn. (1), the diffusion rate cannot be enhanced without increasing the Lorentz force, 
while 3D effects modify the diffusion process itself (compare eqns. (3) and (6)). In this section we 
quantify the effects of magnetic feedback and derive criteria for the validity of the 2D model. 

We assume the stagnation point flow given by eqns. (28) with a = 0. The magnetic field can 
be written using eqns. (32), (34), and (37) as 



B = Boo 



2 



^ ' 27* _ -|N _ ^^^27* 



1 - -^— (e^T* - 1) 



e 



(48) 



47tdo ' ' 4L2 

We evaluate the importance of feedback by following the force on a fluid element over the mixing 
time t^. This probably overestimates the effects of magnetic feedback, because the coherence time 
T of any particular realization of the flow is expected to be less than t^. Thus, the constraints on 
the turbulence which we derive are likely to be conservative. In fact, mixing appears to take place 
in fully self consistent models of Alfven wave turbulence; Maron &: Goldreich (2001). 

In what follows, it is useful to remember that the parameter 2'jtdo, the ratio of the classical 
ambipolar diffusion time to the eddy turnover time, is large; see eqn. (39). We will sometimes use 
the inverse of this quantity as an expansion parameter. 



5.1. Magnetic Pressure Forces 

We estimate the deceleration of an element of fluid by magnetic pressure in time t^. The x 
component of magnetic pressure force Fm is 

"'" = -4l7- <""> 

The deceleration Av^ of a fluid element over a time t is 

1 /"* 

AvP{xo) = - / dsFra{x{s),s), (50) 



P Jo 

where x{s) = xqc"'''* is the position of the fluid element at time s and xq is its original position. 
Substituting eqn. (48) into eqn. (49) and integrating eqn. (50) to t = t-y yields to leading order in 
27tdo 



15 



where we have assumed xq/L ^ 1. 

The average velocity v{xo) of the fluid element over this time is 



1 Z"*^ 
v{xo) = — dsv{x{s),s), (52) 

h Jo 

where x{s) is once again the Lagrangian position of a fluid element. For the stagnation point flow 
(28), 

V = ^, (53) 

to leading order in 2^tdo. Magnetic feedback on the flow is unimportant if Av^ /v < 1. Combining 
eqns. (51) and (53), we derive 

Using eqns. (12) and (39), eqn. ( 54) can be rewritten as 

^ = Ii(^y'\J2^). (55) 

V TA \32TniJ V TdTniJ 



Equation (55) implies an upper limit r^^^ on the eddy turnover time r^^ such that Av /v < 1; 



TA V TA 



-'^^ r\-r\ -{-irx rs ti^rA r\-%T ■|-ttt*vi rw m-v "flTTl f^ t , cti/^K -J-Ko'l- /x/ii-'^^ i 

-2/3 






In' "^ 



\/27 



(56) 



where, to sufficient accuracy, we have replaced t^ by (32r^r„j)"'^' ^ in the logarithmic factor. On this 
timescale the magnetic field is still well frozen to the eddies; from eqns. (12) and (56), 



2\ 1/3 



H..KJ^|(-^) Ki^) 



-2/3 



(57) 



5.2. 3D Effects 

We assume the turbulent motions are in the (x,y) plane, but depend weakly on z; i.e. the 
characteristic wavenumber k along the field is related to the turbulent lengthscale I hy kl <^ 1. 
This quasi-two dimensionality is expected to be a feature of Alfvenic turbulence in strong magnetic 
fields (Strauss 1976, Goldreich & Sridhar 1997). 

Since the field is fairly well frozen in even on scales /, it is very well frozen on the scale k~^, 
and the transverse field B_l is given to a good approximation by 

B^ = 5,(xo,t)^, (58) 
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where xq is the initial position of the fluid element at position x at time t. 

Let us introduce a small parameter e and assume that dz is 0{e) relative to the perpendicular 
derivatives, and that B±/Bz is also 0{e) (this is the so-called reduced MHD ordering; Strauss 1976). 
It can then be shown that the changes in the ambipolar drift terms are of order e^. Therefore, we 
may assume that weak three dimensionality has little effect on ambipolar drift of the vertical field. 

The bent field exerts a tension force which decelerates the fluid by an amount Av'^ . We 
compute Av for a z dependent stagnation point flow model with 

Vx = — ^x cos kz, (59) 

Vy = 7y cos kz, (60) 

a generalization of eqns. (28). The Lagrangian positions are 

x = xoe-^^ ''"'''', (61) 

According to eqns. (58) and (61), the x component of the field is 

Bx = Bz^tkx sin kz. (63) 

The magnetic tension force Fm in the x direction is 

F„ = -1 — ^ = -^-ftk^x cos kz. (64) 

47r oz Ait 

We set z = and follow a procedure similar to the derivation of eqn. (55), integrating F^^ along 

the path of a fluid element from t = to t = t^. We approximate Bz by i?00) which overestimates 

Fm- The result to leading order in (27trfo)~^ is 

Av^=^-f'dtF^=^^^^^. (65) 

Using eqn. (53), the relative deceleration is Av /v 

Av^ k'^vl , , , 

-^ = —fln{2jtdo). 66 

V 27^ 

Equation (66) shows that magnetic tension has little effect on the fluid as long as the Alfven 
frequency along the fieldline is less than the eddy turnover rate by the factor [In (27^,^0)]^ • 

Equation (66) can be used to set an upper limit rl^ax o^ the eddy turnover time such that 
the magnetic field reaches the mixing scale without decelerating the fluid. Proceeding as in the 
derivation of eqn. (56) we find 



r, 



T 
max 



V2 



In (V2kL— 



-1/2 



(67) 



TA kL 

Equation (67) shows that tension forces are less important in long, thin structures, in which kL 
can be much less than one, than they are in flattened structures such as disks. The field is well 
frozen to the eddies as long as TA/Tni 3> 1. 
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6. Application to the Galactic Magnetic Field 

Little is know about interstellar turbulence beyond its gross energetics: the turbulent kinetic 
energy is at or above equipartition with the magnetic energy. Energy injection by a variety of 
mechanisms, combined with nonlinear processes, should lead to a turbulent spectrum over a wide 
range in scales. 

In principle, any weakly ionized interstellar structure which survives for several eddy turnover 
times is a candidate for turbulent ambipolar diffusion. We showed in §5 that the efficiency of turbu- 
lent ambipolar drift can be limited by the back reaction of magnetic forces. In the strictly 2D case, 
the increase in ambipolar diffusion rate is associated with the local buildup of magnetic pressure 
forces, and in the 3D case, by magnetic tension forces as well. We expressed these constraints in 
terms of lower bounds on the strain rates, or inverse eddy turnover times, such that concentration 
of the field occurs before deceleration of the flow. These constraints appear in eqns. (56) and (67). 
Here, we express them numerically. 

The critical parameters are the ratio Tni/rA and the geometrical factor kL. Referring back to 
^2.2 for numerical values, we have 

(68) 

5 X 10"^ cm"^, ^„ = 1.4, Lpc = 
1, S^ = 3, TniJTA = 3.4 X 10~^, and eqn. (56) requires t^/ta < 0.13. If Vt ~ va, turbulence on the 
scale of a tenth of a parsec or less can mix the magnetic field down to the ambipolar diffusion scale. 
On the other hand, eqn. (67) requires t^Jta < 0.09(27r/fcL), which is a more severe constraint, 
especially in a highly fiattened structure. 

There is at least one type of H I structure in which flux freezing appears to be obeyed. Mag- 
netic flelds in H 1 shells are observed to be quite strong, with magnitudes consistent with shock 
compression (Heiles 1989). The same observations suggest that turbulence with Alfvenic or slightly 
subAlfvenic velocities is present. If these shells were not expanding, they would appear to fulflll the 
conditions for fast ambipolar drift, and their strong fields would be counterexamples to the theory. 
However, expansion of the shells at speeds of order 10 - 20 km/s adds new magnetic fiux faster 
than it can diffuse upstream, while on the downstream side the ionization is too high for efficient 
ambipolar drift. Thus, the field in the shells remains large. Recent observations by Heiles (2001b) 
and Heiles & Troland (2001) of cold, moderately dense H 1 regions which have weak magnetic fields, 
Alfvenic random velocities, and no observed association with shells appear to be better candidates 
for fast ambipolar drift. 



Tni 
TA 


= 4.8 X 10-5 ^ ^ 


For ( 


3xample, if rin = 50 cm"'^, rii = 
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7. Summary and Conclusions 

Observations sliow that wiiile interstellar magnetic fieldstrength and gas density are to some 
extent correlated, the fieldstrength is often lower than expected. This suggests that processes 
beyond ideal MHD may play a role. 

The flux to mass ratio is altered by ambipolar drift, but estimates of the ambipolar diffusivity 
v^Tni predict that ambipolar drift is important only in very dense, strongly magnetized gas with 
substantial gradients on small scales. However, it is well established that turbulence can enhance 
the transport rates of quantities such as entropy and angular momentum. This motivated us to 
consider the effect of turbulence on the rate of ambipolar drift. Enhancement by roughly two orders 
of magnitude would explain the observations. 

As a first attempt on the problem, we considered the geometrically restricted case of a straight 
magnetic field, with a transverse gradient, mixed by 2D, perpendicular turbulence. In this situation, 
ambipolar drift is described by a nonlinear diffusion term [eqn. (6)]. For simplicity, however, we 
approximated the diffusivity as linear. 

We showed by a mixing length argument, and then an explicit calculation (§3) that advection 
of the field by a periodic flow reduces its peak value. However, unless the motions are long and 
thin, like streamers, the rate of relaxation is no faster than relaxation by ambipolar drift alone. 
The missing ingredient is stretching and shrinking of scales, which in chaotic flows happens at 
an exponential rate. In §4 we modelled these exponential changes of scale by representing the 
flow as a sequence of randomly oriented hyperbolic stagnation points. With this model, and a 
parabolic initial condition for the magnetic profile, the advection-diffusion equation (7) can be 
solved exactly. The model predicts an exponential increase of the mean diffusion rate with time 
(Figure 1), although with considerable variance from point to point (Figures 2 & 3). The stagnation 
point model predicts that the fleld diffuses on a timescale comparable to the eddy turnover time, 
with only logarithmic dependence on the ambipolar diffusivity itself and on the original gradient 
lengthscale. 

In §5, we estimated the back reaction of magnetic pressure and tension forces on the stagnation 
point flow, including weak three dimensionality. The relative deceleration of the fiuid over one 
mixing time is given for pressure forces by eqn. (54) and for tension forces by eqn. (66). Comparing 
the deceleration time of a fiuid element to the mixing time, we derived upper limits on the eddy 
turnover time such that deceleration is order unity or less within a mixing time [eqns. (56) and 
(67)]. As we showed in §6, these criteria can be satisfied in the interstellar medium without extreme 
assumptions about the size and velocity of the turbulent eddies, although they cannot be wholly 
ignored. These estimates of feedback are conservative in the sense that a fully self consistent 
model of MHD turbulence can still produce fast mixing, as shown by Maron & Goldreich (2001) 
for spreading of a passive scalar. This conclusion may be dependent on geometry, however, as 
concluded by Kim (1997) based on a study of turbulent decay of a coplanar magnetic field caused 
by 2D incompressible motions in weakly ionized fiuid. 
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The outcome of these calculations is that turbulence is likely to have a major effect on the 
magnetic flux to mass ratio in the weakly ionized portions of the interstellar medium, making the 
magnetic field more uniform. The model presented here applies to regions with simple magnetic 
topology, filamentary structure, and no global cross-field flows. It follows that the strength of the 
field is not necessarily a good indicator of the dynamical processes which determine the gas density. 

Although the results presented here are consistent with the conventional wisdom that turbulent 
mixing takes place in an eddy turnover time, it is important to recognize that they are not obtained 
from a full model of turbulence. We plan to extend this work to more realistic models which include 
compressibility, vorticity, and a self consistent treatment of magnetic forces as well as more general 
magnetic geometry. If the stagnation point model holds up in comparison with more complete 
models, it could be useful in other mixing problems. 

Magnetic reconnection can also change the magnetic flux to mass ratio. Lazarian & Vishniac 
(1999) have argued that if the spectrum of interstellar turbulence extends to the resistive scale 
then reconnection takes place at the Alfven speed. We have concentrated here on ambipolar drift 
because it does not require turbulent structure on such small scales; if flux is quickly redistributed 
in the fully ionized portions of the ISM then an alternative process is certainly required. 

This work was initiated during the program on Astrophysical Turbulence at the ITP in Santa 
Barbara in 2000. I am happy to acknowledge useful discussions with Nic Brummell, Dick Crutcher, 
George Field, Carl Heiles, Fabian Heitsch, David Hughes, Eun-jin Kim, Steve Tobias, Tom Troland, 
and especially the referee, Pat Diamond. Material support was provided by NSF Grants AST 
9800616 and AST 0098701 to the University of Colorado and PHY 9407194 to UC Santa Barbara. 
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